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ABSTRACT 

We study the physical properties of 216 z ~ 2.1 LAEs discovered in an ultra-deep narrow-band 
MUSYC image of the ECDF-S. We fit their stacked Spectral Energy Distribution (SED) using Char- 
lot & Bruzual templates. We consider star formation histories parametrized by the e-folding time 
parameter r, allowing for exponentially decreasing (r > 0), exponentially increasing (r < 0), and 
constant star formation rates. We estimated the average flux at 5015 A of our LAE sample, finding a 
non-detection, which translates into negligible Hell line emission at z ~ 2.1. In addition to this, the 
lack of high equivalent width Lya line objects ruled out the hypothesis of a top-heavy IMF in LAEs. 
The typical LAEs of our sample are characterized by best fit parameters and 68% confidence intervals 
of log(M*/M )=8.6[8.4-9.1], E(B-V)=0.22[0.00-0.31], r=-0.02[(-4)-18] Gyr, and age SF =0.018[0.009-3] 
Gyr. Thus, we obtain robust measurements of low stellar mass and dust content, but we cannot 
place meaningful constraints on the age or star formation history of the LAEs. We also calculate the 
instantaneous SFR to be 35[0. 003-170] M yr _1 , with its average over the last 100 Myr before obser- 
vation giving <SFR> 10 o=4[2-30] M Q yr _1 . When we compare the results for the same star formation 
history, typical LAEs at z ~ 2.1 appear dustier and show higher instantaneous SFRs than z ~ 3.1 
LAEs, while the observed stellar masses of the two samples seem consistent. Because the majority 
are low-mass galaxies, our typical LAEs appear to occupy the low-mass end of the distribution of 
star forming galaxies at z ~ 2. We perform SED fitting on several sub-samples selected based on 
photometric properties and find that LAE sub-samples at z ~ 2.1 exhibit heterogeneous properties. 
The typical IRAC-bright, UV-bright and red LAEs have the largest stellar mass and dust reddening. 
The typical UV-faint, IRAC-faint, and high equivalent width LAE sub-samples appear less massive 
(< 10 9 M©) and less dusty, with E(B-V) consistent with zero. 

Subject headings: high redshift galaxies: star-forming galaxies, Lyman Alpha Emitters 



1. INTRODUCTION 

Lyman Alpha Emitting (LAE) galaxies have been ob- 
served at very high redshift, thanks to their strong emis- 
sion lines enabling discovery even without continuum de- 
tection. Many studies have been carried out to investi- 
gate their nature, using photometric and spectroscopic 
data (e.g., |Cowie fc Hu||l998| |Malhotra fc Rhoads||2002[ 



Ouchi et al.||2003l Finkclste in et aLl|2008l |Ouchi et al. 
2008UGronwall et al.||2007[ [Nilsson et al.||2007| |Nilsson 
et al.|2009[ ). A description of the physical nature of LAEs 
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can be o btained from spectral energy distribu tion (SED) 
anal ysis ([G awiser 2009, Walche fet al.| 2011 for reviews 
and Kannappan & Gawiser 2007 for caveats). Multi- 
wavelength surveys allow us to study the SEDs, which 
can be obtained from individual object flux densities, 
when there is enough signal-to-noise, or through photo- 
metric stacking of galaxy samples. At z > 3 LAEs are 
young, nearly dust - free, low-mass objec ts (z 



3.1 Ga- 



wiser et al.||2006a| |Nilsson et al.||2007[ |Gawiser et"aT 



20071 |Lai et al.||20b8MOno et al.||2010aT z ~ 4.5 
stein et al.l|2008i z ~ 5jYuma et al ||2010j z ~ 6 



Finkel- 
Pirzkal 



et al.||2007nOno et al.||2010bMOuchi et aL|[2010| ). Fur 

themore, Gawiser et al. (2007) found that the typical 
LAEs at z ~ 3.1 are characterized by a median starburst 
age of ~ 20 Myr, a low stellar mass of M ~ 10 9 M , 



modest SFR values around 
gible amount of dust (A\ 



< 



' 2 M yr , and a negli- 
0.2), indicating that they 
are galaxies in the early phases of a burst of star forma- 
tion (SF). They also measured the spatial clustering of 
the LAE sample, and inferred that LAEs are hosted by 
the population of dark matter halos with masses greater 
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10-6±8;|. 



From the evolution of 



than log(M min /M Q ) 

dark matter mass w ith redshift (iSheth & Tormen||1999 



Francke et al. 
LAEs can evo 



20081, their results imply that 



ve into Milky Way-type galaxies at z ~ 



(Gawiser et al. 2007). This adds to evidence from SED 
fitting and rest-frame ultraviolet sizes that high-redshift 
LAEs could be building blocks of mor e evolved lower red 



shift galaxie s (Pascarelle et al. 1998[ Ouchi et al. 
Ouchi et al.|[2OT0| : 



2003 



LAEs can therefore occupy the low-mass end of high- 
z galaxy populations, where the bulk of the galaxies lie. 
Other well-studied populations, like Lyman Break Galax- 
ies (LBG) or Sub-Millimeter Galaxies (SMG), represent 
the rare luminous tail of the high-redshift luminosity 
function. However, despite having low mass, LAEs do 
not have particularly high number density, implying that 
they are selected from the low-mass bulk of high-redshift 
galaxies due to special conditions in which the Lya pho- 
tons are produced and are able to escape. These condi- 
tions involve star formation, responsible for the ionizing 
photon production, the spatial and velocity distribution 
of neutral hydrogen in th e interstellar medium ( ISM) and 
the ISM dust geometry ( |Verhamme et al.]|2006| as an ex- 
ample). Large samples are crucial to better constrain 
typical properties of star forming galaxies in general, 
and LAEs in particular. The comparison of LAE prop- 
erties with those of color-selected star forming galaxies 
has shown that Lya emission can also be detected in 
older, more massive, dusty galaxies flShapley et al. 2001 , 
Tapken et"aT]|2007| |Pentericci et al.||2009[ |Kornei eFaT] 



2010). 



Recently Nilsson et al. (2009) claimed that at lower 
redshift [z ~ 2.3), LAEs present more diversity in pho- 
tometric properties and hence a possible evolution from 
z>3toz~2, inl Gyr of Universe life. This is con- 
firmed in thei r recent SED analysis (Nilsson et al. 2011). 



Guaita et al. (2010, hereafter GulO) obtained a sample 
of 250 LAEs at z ~ 2.1 in the Extended Chandra Deep 
Field - South (ECDF-S), one of the fiel ds of MUSYC 
(Multi- w avelength Survey by Yale-Chile, Gawiser et al. 
2006 tf^t- The clustering properties of this sample were 
similar to LAEs at z — 3.1, implying that they also serve 
as building blocks of present-day L* galaxies. The Lya 
luminosity and rest-frame equivalent width (EW) lim- 
its (\og(L Lya ) = 41.8, EWrest-frame = 20 A) adopted in 
GulO were the same as those at z ~ 3.1 studied by Gron- 
wall et al. (2007). This allows us to make an unbiased 
comparison between the photometric properties (GulO) 
and spectral energy distributions (this work) of the two 
LAE sample, investigating the LAE evolution claimed by 
Nilsson et al. (2009, 2011). 

Section ^2] of this paper describes the LAE sample and 
sub-samples based on rest-frame NIR flux density, the 
rest-frame UV slope, represented by the B-R color, rest- 
frame UV magnitude, equivalent width. In Section f|3j we 
describe the stellar population synthesis templates used 
in the SED fitting. We present the results in Section S|4j 
includingthe comparison with SED results at z ~ 3.1. In 
Section |5]we discuss the implications of the results of the 
SED fitting of our sample of z ~ 2.1 LAEs and its sub- 
samples. We summarize the inferred physical properties 

12 http: / / physics.rutgers.edu/~gawiser/MUSYC 



of z=2.1 LAEs in Section [|6] 
We assume a AC DM cosmology co nsistent with 

" their Tab. 2), 
0.26, Ov = 0.74, 



WMAP 5-year results punkley et al.|2009 
adopting the mean parameters it 
H n = 70 km s- 1 



Mpc~ 



= 0. 



2. THE SAMPLE 

2.1. Stacked Multi-wavelength Photometry 

We used the sample of 250 LAEs found in the deep 
narrow-band image of ECDF-S by GulO at z ~ 2.1. The 
ten criteria used in the selection of Lya emitting galaxies 
are listed in §4 of GulO. Morphol ogy analysis of the 250 
z ~ 2.1 LAEs (Bond et al. 2011[ ) revealed the presence 
of 34 extended objects with sizes larger than ~ 0.6" in 
the HST/ACS V broad band (point spread function of 
0.1"). This size corresponds to a half-light radius of 2.5 
kpc at z ~ 2.1. Previous studies on the morphology of 
Lya emitting galaxies hel ped us in decidin g to exclude 
these objects in this work. |Bond et al.|2009 also analyzed 
morphologically a sample of z ~ 3.1 LAEs, using HST 
broad band images. They found that the typical high 
redshift LAE was barely resolved at the ~ 0.1" angular 
resolution of the HST/ACS U-band image and that very 
few were larger than 1.5 kpc. In Bond et al. (2011), some 
extended sizes were found for 1.5 < z < 3 star forming 
galaxies, but objects with half-light radius bigger than 
2 — 3 kpc were very rare at z ~ 2. The median half- 
light radius of z ~ 2.1 LAEs is, instead, 1 kpc and the 
majority of the LAEs are characterized by sizes smaller 
than 2 Kpc. 

These extended objects, if real LAEs, must consti- 
tute a separate population not seen at z ~ 3.1 and are 
worth studying separately. Not yet having spectroscopy 
of them to verify their redshift, we prefer the conserva- 
tive approach of excluding them from SED analysis. The 
majority is characterized by low Lya luminosity and lim- 
ited equivalent width. Excluding them from the sample, 
the photometric results presented in GulO are still valid. 
This is mainly due to the fact that we had reported me- 
dian typical properties and also because number density, 
equivalent width, luminosity distributions and star for- 
mation rate estimations were derived for the half of the 
sample complete in terms of EW. The clustering results 
obtained for the sample of 216 objects is also consistent 
with those published in GulO within the errors. The 
low-EW, low-Lya luminosity part of the sample is surely 
composed by real LAEs, but it is also the most likely 
to be affected by low-redshift interlopers. A careful in- 
spection of ground and space-based data is needed, but 
spectroscopy is the most efficient way to disentangle their 
nature. The median stacking approa ch fr om the observed 
optical to IR we adopt below (see [3.1 for a discussion 
of the use of median statistic) should minimize the in- 
fluence of contaminants unless they are numerous, and 
we indeed observed that the SEDs of the full sample and 
of the sub-samples do not change significantly upon re- 
moving these extended objects. However, we removed 
them from the full sample since we do not think they are 
representative of the typical LAE at z ~ 2.1. Our final 
z ~ 2.1 LAE sample is th us composed of 216 objec ts. 



pro- 



Using the SExtractor (|Bertin &; Arnoute 
gram in dual image mode, we estimated the aperture 
fluxes of the narrow-band detected galaxies in all the 



1996) 
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broad bands (Tab. ffl) of the MUSYC survey sho wn in 
Fig. [T] These inclu ded the U, B, V, R, I, z optical (|CaF 



damone et aL [2010 ) , J,H,K near-infrared ( Taylor et al 
20098, and lfeC~3 .6, 4.5, 5.8, 8.0 infrared bands 



iDariien et al. 2011). The apertures in each band were 
optimized by calculating the curve of growth of bright 
point sources. This led us to choose 1.2" diameter aper- 
tures for all the optical bands, 1.6" for J, 1.8" for H, 
2.0" for K, 1.6" for 3.6 fim, 1.8" for 4.5 ^m, 2.0" for 5.8 
/im, and 2.2" for 8.0 /im. We used the fractional signal 
inside the optimal apertures for point sources to trans- 
form the aperture fluxes to total fluxes. This assumed 
that z ~ 2.1 LAEs were unresolved in ~ 1" seeing, which 
our morphology analysis showed to be an excellent ap- 
proximation. To test the accuracy, we looked for matches 
between our LAE sample and version2 of the SI MPLE 
IRAC catalog in the ECDF-S ( |Damen et al.|2011[ ), com- 
posed by ~ 45000 objects. For the 85 common objects 
(34% of our sample) we could compare the SIMPLE pho- 
tometry with ours. The median fluxes agreed well with 
a maximum difference of 15% in the 8.0 /iin band image, 
which is characterized by the worst resolution. 

We are interested in recovering the typical observed 
SED of z ~ 2.1 LAE population. We stacked the LAE 
measured fluxes (including formal non-detections and 
negative fluxes) by computing the median flux in each 
band. Some objects were excluded from the stacking in 
the NIR and IR bands. The NIR MUSYC images cover 
a slightly smaller area on the sky than the optical ones, 
and in the NIR flux stacks we excluded from the 216 LAE 
sample the galaxies outside this NIR area (3 objects were 
excluded from the flux stacks in J and 47 in H band). 
In the IRAC 3.6 /im image, we selected "clean" regions, 
where the photometry of the objects was not affected 
by neighbors; we also excluded the galaxies with one or 
more than one IRAC detected objects in a searching ra- 
dius greater than 1", but smaller than 4" (~ twice the 
IRAC point spread function), and 3 objects that were not 
properly de-blended. The IRAC fluxes from all 4 bands 
were then obtained from the median among 114 objects. 
The LAEs' own properties played no role in the classi- 
fication of "clean" regions, so the resulting photometry 
should be unbiased. 

The errors on the stacked fluxes were estimated as 
follows. For each data point, there are two sources of 
uncertainty: the photometric error on the measurement 
and the sample variance. The former was determined 
by the quality of the observations and the instrumenta- 
tion, while the latter quantified the error introduced by 
assuming that the properties inferred from one sample 
were a good representation of the general properties of 
LAEs at this redshift in the Universe as a whole. We took 
sample variance into account by means of the bootstrap- 
ping technique (B. Efron, 1979, Bootstrap Methods). 
By Monte-Carlo sampling with replacement our catalog 
of 216 LAEs, we built a large number of catalogs of the 
same size as the original and we computed the sample 
variance error as the scatter in the median fluxes in each 
band. This technique was robust when the number of 
realizations was large enough that the estimation of the 
mean of the median fluxes approached the flux in the 
"real" catalog and when the scatter did not change if 
the number of realizations was increased. Our choice of 
using 100,000 Monte-Carlo realizations amply met both 



these requirements. We ignored the formal aperture flux 
uncertainties while doing the bootstrap, where sample 
variance dominated. The zero-point uncertainties were 
then added in quadrature to the sample variance error. 
The zero-point uncertainty was 5% for all the bands, but 
10% for U and H bands, due to the known difficulties 
in calibrating a blue filter and because the ff-band im- 
age is a combination of smaller fields which do not have 
exactly the same photometric zero-point. The disagree- 
ment between ours and SIMPLE photometry could be 
due to a zero-point offset, so we were taking into ac- 
count that with this additional 5% uncertainty. The ob- 
served (N)IR bands were characterized by higher total 
uncertainty than those observed in the optical bands, 
due primarily to sample variance between the bootstrap 
samples. On average the sample variance contribution to 
the total error budget was ~ 80%. 

2.2. Sub-samples of LAEs 

We are interested in the typical SED of the whole sam- 
ple of LAEs at z ~ 2.1, but also in possibly detecting di- 
verse LAEs populations within the sample. Lacking suf- 
ficient signal-to-noise to fed the SEDs of most individual 
objects, we separated the sample into sub-samples based 
on photometric properties and stac ked their flux es fo llow- 
ing the method described in §2.1| and then in §3.1[ We 
fit the SEDs of the following sub-samples, whose specifics 
are listed in Tab. [2j 

1) The full sample of 216 z ~ 2.1 LAEs (FULLsample 
hereafter). 

2) We split the full sample into two based on the 3.6 
/im flux density. Lai et al. (2008) separated the z ~ 3.1 
LAEs in two sub-samples, based on the IRAC 3.6 /im 2a 
detection (f3.6/^m > 0.3/iJy) and non-detection. They 
observed that the IRAC-detected sub-sample (L08det 
hereafter) was composed by older (160 Myr) and more 
massive (10 10 M ) LAEs than the IRAC-undetected 
(L08und hereafter) objects. Our flux cut at is.efim = 
0.57/iJy corresponds to the same luminosity cut, blue- 
shifted to z ~ 2.1, allowing direct comparison with their 
results. Only objects in "clean" IRAC regions were con- 
sidered for these samples; the resulting IRAC-bright sub- 
sample contains 47 LAEs and the IRAC-faint sub-sample 
contains 67 LAEs. The separation based upon the 3.6 
fim flux was made just among the objects located in the 
IRAC "clean" regions. In the "dirty" regions the 3.6 /im 
flux estimation is not reliable due to effect of neighbor- 
ing sources; this did not affect the fluxes for any of the 
non-IRAC bands, which were accepted for the analysis 
on the full field. 

3) Fig. 6a of GulO plotted the observed R band (rest- 
frame UV) magnitude versus the B-R color (proportional 
to the rest-frame UV slope) and observed a red branch 
at R < 23. We are interested in fitting these galaxies' 
SED to infer their nature. Therefore we identified a sub- 
sample on the basis of B — R color. Approximately 15% 
of the full sample (34 LAEs) is characterized by B — R > 
0.5, red-LAEs hereafter. 

4) We defined two sub-samples of LAEs split at R — 
25.5. The continuum-bright (UV-bright hereafter) sub- 
sample of 118 LAEs enables a direct comparison with the 
SED parameters of R < 25.5 " BX" , star forming ga laxies 
in the same range of redshift ( Steidel et al.| [2004). The 
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remaining 98 LAEs are classified as UV-faint. 

5) We used the estimations of L^ya and EW from 
GufO to define a complete LAE sample of If 9 objects 
with L^ya > fO 42 ' 1 erg sec -1 . We investigate the SED 
dependence on equivalent width, dividing the complete 
sample at the median EW res t_/ rame of 66 A, obtaining 
60 high-EW (EW rest -/rome > 66 A) and 59 low-EW 

(EWrest-frame < 66 A) SUb-SampleS. 

We investigated the possible overlaps among the sub- 
samples in Tab. [3j Just 34% (16/47) of the IRAC- 
bright galaxies are also red LAEs, but 89% (42/47) of 
them are also UV-bright. Among the 118 UV-bright ob- 
jects, 42/118 (36%) are IRAC-bright, 16/118 (13%) are 
IRAC-faint, and the remaining 60/118 (51%) are not in- 
cluded in cither IRAC sample due to nearby neighbors. 
Around 50% of the UV-faint objects also belong to the 
sub-sample of IRAC-faint. Therefore we expect some 
similarities between IRAC-bright and UV-bright and be- 
tween IRAC-faint and UV-faint SED results, but none of 
the sub-samples have sufficient overlap to be effectively 
identical. One exception are the blue-LAEs, which by 
virtue of containing ~ 85% of the total sample also con- 
tain the majority of most sub-samples. We expect re- 
sults from blue-LAEs to be very similar to those for the 
FULLsample. 

2.3. BX galaxy comparison sample 

As a means of comparing the SED properties of star 
forming galaxies at z ~ 2, we fit the SED of a sam- 
ple of spectroscopically confirmed color-selected ("BX" 
hereafter) galaxies found in the MUSYC ECDF-S cata- 
log. We selected these galaxies using a two-color diagram 
including U—V, V — R colors and R band magnitude, 
as in Steidel et al. (2004). The criteria took advantage 
of the decrement in the spectrum due to neutral Hydro- 
gen absorption at wavelengths just bluer than Lya. We 
first generat ed rest-frame s pectra of star forming galax- 
ies using the |Bruzual|2003 code at ages between 2.5 Myr 
and 3 Gyr. W e then added a red dening contribution 
parametrized by |Calzetti et al. 2000 law. We evolved the 
resulting spectra with redshift and calculated the U — V, 
V — R colors at every redshift. The reddening amount 
was the dominant parameter in widering the V — R range 
for star forming galaxies. For our filter curve configura- 
tion, the U — V, V — R space that selected redshifts in 
the range 2.0 < z < 2.7 was defined by 



V-R > -0.15 
U - V > 3.45(U - R) + 0.32 
U - V < 4.46(F - R) + 1.6 
V - R > 0.03(U - V) + 0.24. 



(1) 

(2) 
(3) 
(4) 



We also asked for magnitud es 23 < R < 25.5. By def- 
inition ( Steidel et al.|[2004[ ) BX galaxies are bright in 
the continuum; the R > 23 requirement allowed us to 
exclude galactic stars, which are located in the same 
region of the two-color diagram (GulO Fig. 7b). Just 
like Lyman Break galaxies at z ~ 3, BX galaxies could 
show Lya emission line in their spectra, depending on 
their interstellar medium properties. Among the ~ 84000 
sources in the MUSYC catalog, we considered the ~ 6000 



BX galaxies present in our narrow-band detected cata- 
log. About 200 of them were characterized by spectro- 
scopic redshift. We then chose the BXs charac terized by 
2.0 < z spec < 2.2 from the GOODS Mastercat PPj and 
MUSYC databases. We excluded galaxies that showed 
an IRAC 3.6 (im detected source within an annulus of 
1" internal and 4" external radius from their position. 
As a result, this sample is characterized by 16 objects 
with comparable photometry to our LAE sample. Due 
to Lya selection, most LAEs are faint in the continuum. 
The UV-bright LAE sub-sample is expected to overlap 
with BX galaxies in this property and so we investigate 
if it is reflected in the SED properties. 

3. SED FITTING METHODOLOGY 

We used the S. Chariot & G. Bruzual (2010, private 
communication, hereafter CB10) code to compute the 
stellar population synthesis models. The updated ver- 
sion of the original code includes the treatment of the 
thermally pulsating a symptotic giant br anch (TP-AGB) 
stars (prescription of Marigo & Girardi 2007). We used 
solar metallicity, for which the stellar population model 
was described as "very good" in the Bruzual & Char- 
lot (2003) documentation. In addition to that, the con- 
tinuum model SED did not show significant differences 
in the shape when changing the metallicity value. In 
order to be able to compare with previous works we 
adopt ed a Salpeter initial mass function (IMF, 
1955[ ), distributing stars from 0.1 to 100 M , 
$(M ) oc M" 2 - 35 . However, a Chabrier IMF (|Chabrier 



Salpeter 
: bllowing 



2003[ ), for example, being characterized by a flatter slope 
might be more appropriate for deriving an absolute value 
of the stellar mass or SFR, while a Salpet er one can over- 
estim ate these quantities by a factor 1.8 (Papovich et al. 
20111 



This choice of population synthesis model permitted us 
to add a star for mation history (SFHs) parametrized b y 
the r parameter ( Maraston et al.||2010| Lee et aL| 2010), 

SFR(t) = ^e-*. (5) 

N was the normalization of the model which gave the 
observed stellar mass as a function of the star formation 
history. Equation |5]) allowed the star formation to be ex- 
ponentially declining (positive value of r) , exponentially 
increasing ( negative value of r) or constant (r — > oo). As 



described i n [Maraston et ah (2010), an exponentially in- 
creasing SFR could be a good representation of the star 
formation history of a starburst galaxy at z ~ 2, the 
epoch of peak cosmic star formation density. 
For |t| >> t, the SFH approaches the constant star for- 
mation rate, 

N „-± _ N_ 

FT' 

For < r << t, the model describes a single burst, with 
SFR — > at late times, 
limt =0 

For t < — t (r > t), equation ^ can be approximated as 
a linear increase (decrease) in time, e.g., lim±_^. j^ e~ = 

13 http:/ / archive. eso.org/cms/eso-data/data-packages/goods- 
fors2-final-data-release- version-3 . 

14 http://archive.eso.org/cms/eso-data/data-packages/goods- 
vimos-spectroscopy-data-release- version- 1.0/ index 
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Therefore, this model is able to produce approximately 
constant, linearly increasing/decreasing, exponentially 
increasing/decreasing and single burst star formation his- 
tories. We will adopt the exponentially increasing, de- 
creasing and constant star formation (CSF) as the three 
main scenarios of star formation histories. 

The principal parameters used for the SED fitting were 
the observed stellar mass, M*, the extinction, E(B-V), 
the e-folding time, r, and the age since the beginning 
of star formation, agegp. These were the input param- 
eters that went into the computation of each model's 
SED, to be compared with the data. We defined a four- 
dimensional grid of points, the "parameter space", by 
sampling several tens of values for each parameter. We 
used a logarithmic spacing for agesF> which varied be- 
tween 0.005 Gyr and 3 Gyr, also for t, which varied be- 
tween 0.005 Gyr and 4 Gyr, as well as for the observed 
stellar mass, whose range was 10 6 — 10 12 M Q . The age 
of the Universe at z ~ 2.1 is ~ 3 Gyr. For E(B-V), we 
adopted a linear spacing between 0.0 and 0.6. Tab. [4] 
summarizes the parameter grid. 

We took the luminosity densities (erg s _1 A -1 per unit 
of solar luminosity) output by CB10, applied the Calzetti 
dust extinction law and a correction for the intergalac- 
tic neutral H ydrogen absorption using the Madau law 
( Madau 1995 ). We then convolved the flux densities with 
the filter transmission curves (Fig. [IJ. 

For each model (identified by a point in the grid) we 
computed a x 2 value as described below. To find the best 
fit model and compute confidence levels of the fit param- 
eters, we studied how the probability density changed 
along the grid. This was done by computing, 



13 

E 
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(6) 

where i ran over the 13 observed bands, /™° del and f°\ s 
were respectively the theoretically predicted (for each 
model) and observed flux densities, and (Tj was the un- 
certainty on the i-th data point. 

Once the \ 2 was computed for each combination of 
parameters in the grid, we identified all models that had 
X 2 within the 68% confidence level from the expected x 2 
for 9 degrees of freedom. We rescaled the delta-x 2 corre- 
sponding to the 68% confidence level, based on the x 2 of 
the best fit. This corresponded to the best fit in a Fre- 
quentist approach. For each parameter, we reported the 
68% confidence level range as the parameter uncertainty. 

CB10 produces the galaxy continuum luminosity den- 
sity. The spectral energy distribution of a galaxy can also 
be influenced by nebular continuum emission, which can 
make a significant difference when the galaxy is younger 
than 5 Myr, and nebular emission lines, which can affect 
the results for ages up to ~ 50 Myr. Hot stars in active 
star-formation regions provide a large number of ionizing 
photons. They excite neutral Hydrogen (and also Helium 
and Oxygen) producing nebular emission lines, but also 
a gas continuum luminosity that is proportional to the 
Lyman continuum {X re st— frame < 912 A) photons. The 
relevant nebular emission lines produced in a star form- 
ing galaxy at z — 2.1 are the Lya emission itself which 



enters the observed U band (we are taking into account 
its contribution, by subtracting Lya flux from the U con- 
tinuum as in GulO), [Oil] A3727 in the J band, H/3 and 
[OIII] A4959, A5007 in the H band, Ha in the K band, 
and Paa in the 5.7 /im band. 

The inclusion of the nebular emission is beyond the 
scope of this paper. Here we adopted a standard IMF 
(Salpeter, 1955), a reference metallicity, Zq, and a pure 
continuum in the SED model to be able to explore a more 
complicated SFH. We will include a proper tr eatment of 
the nebula r emission in a subsequent paper (Acquaviva 
et al.|2011>. In any case, the nebular continuum is wave 
length dependent and affects mainly the rest-frame opti- 
cal. As its effect falls within our error bars, we can ignore 
its contribution to the SED at this stage. In addition to 
the nebular continuum, at z ~ 2.1 the nebular emission 
lines also enter observed bands that are characterized by 
high uncertainties (see Fig. [2]). Therefore, they are un- 
likely to affect the SED results we show in §4. We checked 
that when including nebular continuum and emission 
lines in a constant star formation model, the resulting 
SED fitting parameters were consistent within the 68% 
confidence level for the FULLsample, analyzed without 
nebular emission. In future analyses, better NIR pho- 
tometry is needed to better constrain the age of LAEs. 
At that point the nebular emission inclusion would be 
important, because models with only stellar continuum 
can appear older than those with the proper nebular gas 
treatment. 

Nebular emission lines can also come from Helium. 
Hell emission at A = 1640 A enters the observed B 
band. Since Hell A 1640 emission is weaker than Lya 
(it is expected to have EW of just few A, Schaerer & 
Vacca 1998, Brinchmann et al. 2008) and Lya is a small 
part of the U band flux, we can be confident that the B 
band flux is not significantly affected by Hell A1640. We 
searched for Hell line detection {Xobserved- frame — 5028 
A at z — 2.066), estimating the flux density of our LAEs 
in the NB5015 image of the ECDF-S (Ciardullo et al. 
2010 in prep). The stacked flux density of the entire sam- 
ple revealed non-detection of Hell. The average NB5015 
flux is -0.0075 /LtJy with a scatter of 0.058 /iJy. This can 
be translated into a mean EW(HeII A1640) of -1.7 ± 1.0 
A, where the errorbar gives the standard deviation of 
the mean derived from the observed scatter. A small 
equivalent width may be expected for young starburst 
galaxies. EW(HeII A1640) would be larger than 5 A in 
the case of an IMF characterized by a higher number of 
massive stars than a normal one, such as in the case of 
a top-heavy IMF. This unusual IM F would also produce 
objec ts with EW(Lya) > 240 A ( iMalhotra fc Rhoads 
2002), which we did not find in significant number in 
our sample. We neither found objects with EW(HeII 
A1640) > 5A at 95% confidence level. In the distribution 
of EW(HeII A1640) of our z ~ 2.1 LAEs there is an equal 
number of objects characterized by strong positive and 
negative EW(HeII A1640), implying that all of these ob- 
jects could be affected by photometric errors rather than 
true strong Hell emissions. Even if only spectroscopy 
can confirm this, we do not see evidence either for weak 
Hell emission from the population or for a significant 
number of individual objects with strong Hell emission. 
As a result, the non-detection of Hell flux density argues 
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against the hypothesis of a top heavy IMF in the LAE 
population. 

We reported the input parameters of the SED, M* 
in M , E(B-V) in magnitudes, r, and agesF in Gyr, 
as free parameters of the SED fit. We also computed, 
as derived parameters, the mean stellar population age, 
t*, in Gyr, the instantaneous SFR, the averaged SFR 
over the last 100 Myr of galaxy evolution, <SFR>ioo, 
the SFR corr (UV), obtained from the dust-corrected rest- 
frame UV flux density, in M Q yr" 1 , and the agesv/T 
ratio. The t* parameter was determined to be the aver- 
age of the age of the stars in a period corresponding to 
the star formation phase, 



t* = age SF 



J a9£sF t x SFR{t)dt 



(7) 



The instantaneous SFR was calculated using Equation[5j 
while the averaged SFR was defined as 



< SFR > 100 = 



lOOMyrs 



SFR(t)dt 



lOOMyr 



(8) 



We chose 100 Myr, because this was the minimum SF 
time for which t he SFR corr (UV) estimator was robust 
( Kennicutt|1998 ). We calculated the latter from the rest- 
frame UV (1500-2800 A, observed R band) luminosity 
density (Equation 1 in Kennicutt, 1998) and assuming 
the E(B-V) range allowed by the SED fit for the dust- 
correction (GulO). Kennicutt's formula assumed a con- 
stant star formation rate and that the galaxy was at least 
100 Myr old. For a 10 Myr old galaxy, the inferred SFR 
was estimated to be roughly half the real one. Therefore, 
for a CSF model, <SFR>i o was directly comparable 
with SFR corr (UV). However if the galaxy was younger 
than 100 Myr, <SFR>i 00 and SFR corr (UV) were under- 
estimations of the on-going SFR. Also, for an exponen- 
tially decreasing SFH, we would already expect that the 
instantaneous SFR is lower than <SFR>ioo- Instead, 
assuming an exponentially increasing SFH, we would ex- 
pect an instantaneous SFR larger than <SFR>i o- 

The agesF /t derived parameter was determined to in- 
dicate the number of e-foldings that the SFH has un- 
dergone. In the case of exponentially decreasing SFR, a 
large value of agesF/f implied passive evolution with a 
negligible ongoing SF an d that the majority o f the stars 
in the galaxy were old (Grazian et al. 2007). A large 
value of agesF /r in the exponentially increasing SFR 
case meant the agesF was a poorly constrained parameter 
since the implied old stellar population was a very small 
fraction of the mass. The total "stellar" mass of a galaxy, 
which includes the mass of stellar remnants and outflows 
but does not account for return of outflows into future 
generations of stars, can be calculated analytically from 
the integral of SFR(t), Mgal(age SF )=/ O!?esF SFR(t) dt. 
This turned out to be 10-20% higher than the observed 
stellar mass, M*. 

3.1. Median statistics in SED fitting analysis 

As most individual LAEs in our sample had insuffi- 
cient S/N to allow meaningful SED parameter estima- 
tion, the only way to learn about the majority of the 
sample was to analyze their stacked SED. A few galaxies 



with sufficiently bright continuum could be fitted indi- 
vidually, but, doing this would bias the results in fa- 
vor of the properties of objects with bright continuum. 
The stacking approach was also adopted by Lai et al. 
(2008) and Gawiser et al. (2007) to fit the sample of 
ECDF-S LAEs at z ~ 3.1. Outliers in the LAE dis- 
tribution and other contaminants are characterized by 
observed SEDs different from the typical one. As the 
mean stack of the individual SEDs would be sensitive to 
outliers, we stacked the measured LAE fluxes (including 
formal non-detections and negative fluxes) by computing 
the median flux in each band. The analysis of the full 
sample reveals the property of the typical LAE galaxy; 
to investigate heterogeneity in physical properties that 
might be erased through our stack procedure w e sta cked 
sub-samples based on photometric properties (S2.2). 

We created 200 model galaxy SEDs with M77"^E(B- 
V), r, and agesF parameters normally distributed around 
M* = 2.9 x 10 9 M Q , E(B-V) = 0.21, r = 0.1 Gyr, and 
agesp = 0.053 Gyr. These represent typical z ~ 2.1 
LAEs in our observed sample. We added 40 (15%) simu- 
lated z ~ 2.1 galaxies with very different properties with 
respect to the typical ones, and 20 (8%) low-redshift in- 
terlopers. We generated a median stack of SED fluxes 
and associated error bars to the fluxes in the same way 
used for the real data, including the sample variance er- 
ror, evaluated via the bootstrapping technique, and the 
photometric errors. The results of the SED fit were 
M* = 2.82[1.78 - 3.16] x 10 9 M Q ; E(B-V) = 0.22[0.20 
- 0.28]; ; r = 0.087[(-2.952)- (2.759)] Gyr; age SF = 0.057 
[0.023 - 0.143] Gyr. 

We also generated 200 model galaxy SEDs start- 
ing with input parameters uniformly distributed in the 
ranges 8.3 < log(M,/M ) < 9.9, 0.00 < E(B-V) < 0.29, 
-2.3 < log(r) < 0.6, -2.05 < log(age SF ) < -0.64. The 
median input parameters were log(M*/MQ) = 9.1, E(B- 
V) = 0.16, log(r) = -0.87, log(age SF ) = -1-34. Then we 
included outliers as objects characterized by much lower 
and much higher stellar mass at z ~ 2.1 and located 
at lower redshift. We stacked the corresponding fluxes 
calculating their median and fitted the resulting model 
SED. We recovered the four parameter ranges allowed 
within the 68% confidence level, log(M*/M ) = [8.95 - 
9.45]; E(B-V) = [0.08 - 0.25]; log(r) = [(-2.3)- (0.6)]; 
log(agesF) = [(-1-9) - (-0.79)]. They correspond to the 
68% of the input ranges. The median stacked SEDs are 
able to reproduce the typical properties of a sample even 
in the presence of outliers and low-redshift interlopers. 

In addition to that, to test the reliability of our code 
results, we ran it on Mock catalogs, generated from CB10 
models fixing particular values of the free parameters. 
In all the tested cases (r > 0, r < 0, agesF > \t\, and 
agesF < |t|) we could recover the original parameters of 
the models within the 68% confidence level. 

4. RESULTS 

We fit the observed stacked SEDs of the full sample of 
z ~ 2.1 LAEs and those of the sub-samples enumerated 
in 



4.1. FULL sample 

In Fig. [2] we show the observed SED of the full sam- 
ple of z ~ 2.1 LAEs. The best fit model SED is ob- 
tained for an exponentially increasing star formation 
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model. It is characterized by x 2 = 12.3 for 13 data 
points for the combination of parameters labelled in red 
in the figure. The exponentially decreasing and constant 
star formation histories (SFH) produce almost equally 
good fits. Four models at the edge of the 68% confi- 
dence level (x 2 —21.94) are also shown. They are cho- 
sen to be characterized by the highest allowed M* (~40 9 
M ) with E(B-V)=0.07 or E(B-V)=0.1 and lowest al- 
lowed M, (2.82 x 10 8 M ) values with E(B-V)=0.16 or 
E(B-V)=0.17. The other models allowed at the 68% con- 
fidence level have SED shapes within these extremes. 

Tab. [5] shows the best fit and 68% confidence level 
results for the entire parameter space and for the three 
star formation histories separately. In addition to the 
free parameters, it also reports the instantaneous SFR 
and <SFR>ioo- For the best fit model, we obtain 
an instantaneous SFR=35[0. 003-170] M yr _1 and a 
<SFR> 100 =4[2-30] M yr" 1 . 

M* and E(B-V) are quite well constrained while there 
is a wider uncertainty in agegF and essentially all values 
of r are allowed at 68% confidence. This is more easily 
seen in Fig. [3j where we show the contours of the allowed 
parameter space regions for the three main SFHs. In the 
figure we can see the trends of E(B-V) versus agesF (top 
left panel) and versus r (top right panel). Reddening and 
agesp are mildly degenerate; the shape of the SED can 
be similar if the reddening is higher or if the spectrum 
is older. The middle right panel shows M* versus E(B- 
V), while the bottom right panel shows r versus agesF- 
For the exponentially decreasing SFH, the best fit agesF 
and t are both close to 10 Myr. However, the allowed 
parameter space is uniformly covered in the agegF and 
r directions. For the exponentially increasing SFH, the 
best fit agesF and |r| are both close to 20 Myr. High val- 
ues of \t\ together with high agesF values are disfavored 
at the 68% confidence level. Models with agesF up to 1 
Gyr are allowed for SFHs corresponding to |r| below 100 
Myr. Models with agesF up to 3 Gyr are allowed for |t| 
around 100 Myr. 

4.1.1. Comparison with z ~ 3.1 LAEs 

In Tab. |6]we compare the SED fitting parameters with 
others p ublished for LAEs at r edshifts 2 and 3. Pre vious 



work b y [Nils son et al. | ( |2007| , hereafter N07 and jOno, 
et al. ( |201Uap , hereafter U10 assumed a constant SFR 



history, so we focus on that model for comparison. In 
the case of N07, we report the average behavior of the 
z ~ 3.1 fitted galaxies; in the case of O10, we report the 
parameter range of their AT-detected galaxies and the 
results from the st ack of the AT-undetected sources. IGa-l 
wiser et al.] ( |2007[ ), hereafter Ga07, fit the stacked SED 
of the sources with non-detection at 2er in the 3.6 /jm 
IRAC band image. They used a two-population model 
with each population having an exponentially decreasing 
SFH. Nilsson et al. (2009) reported that the prelimi- 
nary SED fitting of their z ~ 2.3 LAEs shows they are 
older and dustier than their z ~ 3.1 sample, noting a 
wider range of possible fitti ng models. Such a result was 



recently expanded upon in [Nilsson et al 
after N10). They estimate the 



(2011, here- 
age of the old and of the 



young stellar population and the mass fraction in the 
young component. As for their entire sample the mass 
fraction in the young component is 56%, we reported 



the young age in the table. While a direct compari- 
son of ages and star formation histories with our work 
is not possible since they use two Simple Stellar Popu- 
lation (single-burst) models and two stellar populations, 
they confirmed the earlier finding of a higher amount of 
dust at lower redshift, and report high stellar masses of 
[5-10] xlO 9 M . We see that our LAEs at z - 2.1 are 
dustier than those at z ~ 3.1 and than the A'-detected 
galaxies from O10, while the observed stellar masses are 
consistent with N07 and Ga07. They are significantly 
higher than the masses of the sample of A'-undetected 
sources from O10 and lower than the range of values ob- 
tained by the same authors for the AT-detected sources. 
The agesF of our sample is consistent with the results 
for A"-undetected sources from Ono et al. (2010), but 
a wide range of values is allowed. Our estimated SFR 
is higher due to the higher inferred dust extinction and 
lies in the range of values found by O10 for their K- 
detected sources. We will compare more directly the Lai 
et al. (2008) res ults w ith our IRAC-bright and IRAC- 
faint samples in §4.2. 1| 

4.2. IRAC bright and faint sub-samples 

The first sub-samples we fit are those selected based 
upon the IRAC 3.6 ^m flux (Tab. In Fig. E] (upper 
left panel) we present the observea and best nt model 
SED for the IRAC-bright and IRAC-faint galaxies. The 
IRAC-faint sub-sample is characterized by a flat rest- 
frame UV-slope and thus by a low dust amount; the ob- 
served (N)IR bands present a higher flux for the IRAC- 
bright sources, implying an order of magnitude higher 
observed stellar mass. This is consistent with Fig. I5J that 
shows stellar mass versus reddening contours of the sub- 
samples at 68% confidence level. These contours do not 
overlap for IRAC-bright and IRAC-faint galaxies. There- 
fore, the IRAC-faint galaxies appear to be significantly 
less massive and less dusty than the IRAC-bright. In 
Fig. [6] the contours indicate the E(B-V) versus agesF 
allowed regions for the sub-samples. While the E(B-V) 
68% confidence level regions show the IRAC-bright typ- 
ical object is significantly dustier than IRAC-faint, the 
agesF contours overlap. 

The IRAC-faint t versus agesF contours (Fig. j\ oc- 
cupy almost the entire parameter space. The fRAC- 
bright contours show a different behavior in the case of 
exponentially decreasing or increasing SFHs. When fit- 
ting the data with an exponentially decreasing model, 
agesF is better constrained (< 20 Myr), while with an 
exponentially increasing model higher agesF values are 
allowed, but just for |r| shorter than 100 Myr. 

Tab. [[[presents results for free and derived parameters 
with 68% confidence levels for the sub-samples. We can 
see that the IRAC-bright galaxies' SED is consistent with 
a < 160 Myr old population, with an <SFR>i o of [10- 
120] M yr -1 and an instantaneous SFR of [40-1310] 
M yr _1 . The IRAC-faint sub-sample is characterized 
by a stellar population of agesF < 1.5 Gyr, with a lower 
<SFR> 100 of [1-15] M yr- 1 . 

4.2.1. Comparison with IRAC detected and undetected 
z ~ 3.1 LAEs 

In Tab. [6] we also compare the IRAC-faint and IRAC- 
bright sub-sample results with the Lai et al. (2008) 
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L08det and L08und results. Lai et al. used the Bruzual 
& Chariot (2003) version of the stellar population synthe- 
sis model in their analysis and assumed a CSF. Our sub- 
sample luminosity threshold matches theirs, allowing for 
a fair test of redshift evolution. We can see that the E(B- 
V) estimation obtained for our sample implies a higher 
dust amount at z ~ 2.1 than at z ~ 3.1. The observed 
stellar masses are consistent between the z ~ 2.1 IRAC- 
faint and z ~ 3.1 L08und sub-samples, but are signifi- 
cantly smaller for the IRAC-bright sub-sample than they 
were for the L08det. Lai et al. (2008) found that L08det 
galaxies were older and more massive than L08und ones, 
but that both categories of galaxies were dust-free. Here 
we find that IRAC-bright are dustier and more massive 
than IRAC-faint, with IRAC-bright requiring non-zero 
dust. 

4.3. Red LAE sub-sample 

In the upper right panel of Fig. Hlwe show the observed 
and best fit model SEDs for the sub-sample of red-LAEs. 
They are significantly more massive and dustier than the 
typical galaxies, as is also seen in Fig. [5j It is interest- 
ing to note that in spite of being dusty and massive, the 
stellar population of the red sub-sample can be as young 
as a few tens of Myr, implying that the population dom- 
inating the SED is still a young one, with a very high 
instantaneous SFR ([1-2880] M yr" 1 , within the entire 
parameter space). In Fig. [6^ we can see that the E(B-V) 
allowed ranges are significantly different for the red with 
respect to typical LAEs, while the agegF ranges overlap. 
Fig. [7] shows that only a few locations in the r versus 
agegF parameter space are disfavored for the red-LAEs. 
However, exponentially increasing models with agesF/|7"| 
ratio close to 1 are ruled out at the 68% confidence level, 
for agesF > 0.1 Gyr. Their SF phase could be as long as 
2 Gyr (Tab. ^. 

4.4. Rest-frame UV bright and faint sub-samples 

The rest- frame UV-bright and faint sub-samples (Tab. 
[2]) are analyzed here. In the lower left panel of Fig. H 
we show that the UV-faint galaxies can be characterized^ 
by a low-dust SED, with a quite flat rest-frame UV. The 
contour plots in Fig. [5]and Fig. |6]also show that the UV- 
bright typical galaxy is more massive and dustier than 
the UV-faint. 

The UV-bright objects are characterized by R < 25.5. 
As with the IRAC-bright sub-sample, brighter contin- 
uum enables better age determination. In fact for the 
typical UV-bright LAE the allowed region of the agesF 
is as narrow as [0.007-0.404] Gyr within the entire param- 
eter space (Tab. ^ and [0.007-0.064] Gyr in the case of 
the exponentially decreasing model (Tab. [8]). However, 
for the UV-faint SED basically the entire agesF range 
is allowed. In Fig. [7] the |r| versus agesF parameter 
space seems to follow the pattern seen for the IRAC- 
bright sub-sample. The exponentially increasing models 
with |t| > 100 Myr are ruled out for agesF > 20 Myr. 

The instantaneous and averaged SFRs can be much 
higher for the UV-bright than the UV-faint sub-sample. 
This is expected because the SFR is proportional to the 
rest-frame UV flux, traced by the observed R band, mul- 
tiplied by a correction for dust extinction. 

4.5. High-EW and low-EW sub-samples 



The high-EW and low-EW sub-sample (Tab. [2) results 
are presented here. The lower right panel of Fig. 4]shows 
their observed and best fit model SEDs. 

In Fig. [5j [6] and [7] we can see that the 68% confidence 
level contours of these sub-samples overlap. However, 
the typical low-EW galaxy shows slightly higher stellar 
masses (10 8 5 -10 9 - 2 M ) than the high-EW (10 7 5 -10 8 - 9 
M Q ). The high-EW LAEs appear to be the least massive 
objects (Fig. [5]). This could be caused by these galaxies 
experiencing one of their first phases of star formation 
and thus not yet having formed as much stellar mass as 
typical LAEs. The E(B-V) allowed ranges seem to be 
consistent for both sub-samples (see also Tab. [7]and[8|, 
revealing that high-EW LAEs can be characterized by a 
moderate amount of dust, even if models implying the 
absence of dust also produce acceptable fits. 

<SFR> 10 o can be five times higher for the low-EW 
than for the high-EW sub-sample. The stellar population 
can be as old as 1.5 Gy for the high-EW, but appears 
younger than 900 Myr for the low-EW galaxy. 

We also examine the dependence upon Lyct luminos- 
ity, splitting the galaxies with LL ya > 10 42 1 erg sec -1 
in half. We do not find significant differences in the 
allowed ranges of the SED parameters. We can com- 
pare our reliable estimates of M* and E(B-V) for the 
Ll^q > 10 42 1 erg sec -1 objects with the recent results 
by N10 at z ~ 2.3 (which reached a factor of 1.8 less 
deep Lya luminosity limit). Looking at the entire pa- 
rameter space we find agreement with their results for 
dust reddening, but our SED analysis implies a stellar 
mass roughly one order of magnitude lower. 

4.6. BX star forming galaxies 

In Fig. [8] we show the observed and best fit SEDs ob- 
tained from stacking a sam ple of ECDF-S star forming 
galaxies at z ~ 2.1 (< |2.3[ ). Within the entire parame- 
ter space, the SED best fit is obtained for an exponen- 
tially increasing SFH (brown solid line). The best fit 
and allowed parameter ranges are log(M*/M Q )=10.0[9.5- 
10.1], E(B-V)=0.1[0.0-0.4], and ages F =0.09[0.01-0.20] 
Gyr. Our data for BX cannot discriminate between dif- 
ferent SFHs. The dotted cyan and solid blue lines overlap 
and represent the exponentially decreasing and CSF best 
fit models. The magenta dashed line shows the Erb et al. 
(2006) SED, assuming the median best fit parameters of 
their sample to build the model spectral energy distribu- 
tion. This is a very good fit for the rest-frame UV part 
of the spectrum, but corresponds to a much higher IR 
flux density, producing a^ 2 ~ 80. 

As a comparison we show the best-fit model SEDs 
(green solid line) for the UV-bright LAEs, selected to 
have the same observed R band magnitude limit as the 
BXs. They present a much flatter rest-frame UV con- 
tinuum, caused by a much lower reddening and/or also 
to a younger stellar population. Our limit in the NIR 
does not allow us to constrain the agesF parameter or to 
discriminate between the two possibilities. The typical 
UV-bright LAE (log(M*/M )=9.1[8.9-9.4] from Tab. [7! 
appears less massive than either our BX or Erb et al. 
(log(M,/M0)=lO.3[lO.2-lO.4]) star forming galaxies. 

In Tab. [Hjwe present the best fit and allowed parameter 
ranges for a CSF model to compare with z ~ 2.1 and 
z ~ 3.1 LAE literature. Although the small BX sample 
is characterized by large sample variance, the typical BX 
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galaxy appears dustier and more massive than our LAEs 
at z — 2.1 and previously studied LAEs at z ~ 3. Their 
higher stellar mass implies a more evolved stage than 
LAEs, perhaps as a consequence of previous phases of 
star formation. LAEs may occupy the low mass tail of 
the BX mass function. 

5. DISCUSSION 
5.1. SED parameter constraints 

Due to the smaller number of objects in the sub- 
samples, their observed SEDs are characterized by higher 
uncertainties than the FULLsample. Although each sub- 
sample contains a significant fraction of objects not in- 
cluded in any other sub-sample, Tab. [3] shows suffi- 
cient overlap that similarities in the SED properties could 
be expected for the typical UV-faint(bright) and IRAC- 
faint(bright) sub-samples. Since we are considering the 
median stacked SEDs, it is worthwhile to study the SED 
fitting parameters of both of them. 

We presented the SED stacking and fitting results in 
Fig. [2] for the FULLsample and in Fig. [4] for the 
four pairs of sub-samples, with the FULLsample best 
fit model SED drawn for comparison. We do not pos- 
sess enough signal-to- noise to fit all, or even most, galax- 
ies individually. However, our sub-sample analysis allow 
us to investigate the heterogeneity of our sample, i.e., 
whether there are differences between the typical sub- 
sample member and the typical z ~ 2.1 LAE. In Fig. [4] 
the observed SEDs already seemed significantly different 
between pairs of sub-samples. In all sub-sample SEDs, 
the rest-frame UV fluxes are characterized by smaller er- 
ror bars, while the N(IR) bands have larger uncertainties 
(Tab. nh. For this reason the physical properties we can 
obtain from the rest-frame UV spectral range are well 
constrained. The solid lines in the figure corresponded to 
the SEDs built using the best fit M*, E(B-V), r, age S F 
parameters listed in Tab. [7] The choice of parameter 
pairs in Fig. [5l [6j and [7] was made to show the con- 
tours of two well constrained, one constrained and one 
non constrained, and two non constrained parameters, 
respectively. 

5.1.1. Stellar mass, M*, and reddening, E(B-V) 

M* and E(B-V) are well constrained. In general, we do 
not see a strong correlation between E(B-V) and stellar 
mass. Looking at horizontal pairs of panels in Fig 
we conclude that the typical IRAC-bright LAE is signi 
icantly more massive and dustier than the IRAC- faint, 
that the UV-bright LAEs are slightly more massive and 
significantly dustier than the UV-faint, and that stellar 
mass and reddening contours overlap for the typical high- 
EW and low-EW galaxies. Looking at the figure along 
the vertical direction, it appears that the red-LAEs as- 
sume the highest stellar mass and reddening values, while 
the high-EW galaxies have the lowest stellar mass and 
dust amount among all the sub-samples. 

5.1.2. Star formation age and e-folding time 

The age since star formation began, agesF 7 and the e- 
folding time, r, are very difficult to constrain with our 
data. One of the reasons why the agesF parameter is 
difficult to constrain is the high uncertainty seen in the 
observed-frame NIR bands, where the Balmer break falls. 



The age of a young stellar population that dominates 
the rest-frame UV continuum (and, presumably, the Lya 
photons) is, thus, represented by the mean stellar age pa- 
rameter, t*, with the overall age since the beginning of 
star formation being very difficult to determine due to 
the low luminosity of evolved stellar populations. When 
considering stacks of observed SEDs with higher S/N in 
the continuum, agesF appears better determined. This 
is, for example, the case of th e sta ck of BX continuum se- 
lected star forming galaxies (j |4.6| agesF = [10-200] Myr). 
In general we see that agegF appears better constrained 
in the exponentially decreasing SFH and in particular in 
the CSF case, due to the loss of one degree of freedom 
via the constraint of M* =age x SFR. When assuming 
an exponentially increasing SFH, the large amount of 
recently formed stars, which dominate the model SED, 
make it nearly impossible to determine the agegp in the 
galaxy. 

In many cases, similarly good best fits can be obtained 
from r > 0, t < 0, or constant SFR models and slightly 
different combinations of M*, E(B-V) and agesF- This 
implies that the S/N in the data does not allow us to 
single out one of them. 

Although our photometry does not allow us to con- 
strain agesF and r, Fig. [7] shows that there are regions 
of the parameter space which are disfavored at the 68% 
confidence level. For the IRAC-bright and UV-bright 
sub-samples, models older than > 20 Myr and with a 
|r| parameter higher than 100 Myr are ruled out. As a 
consequence, t* is required to be < 160 Myr for these 
sub-samples. 

Previous work fitting SEDs of star forming galaxies 
with an exponentially declining SFR also found the star 
formation e-folding pa r ameter difficult to constrain (e.g. 



Papovich et al.1[200 1[) 
set of z ~ '6 LBGs, 



Pentericci et al 
some ot which show 



. (2007) fitted a 
"Lya in emission. 
The agesF and r parameters they found were allowed al- 
most in their entire range, but they were able to constrain 
the agesF over r ratio which they called the evolutionary 
stage of star formation. An agesF/T > 4 was charac- 
teristic of an evolved galaxy, assuming a decreasing SFR 
model (Grazian et al. 2007[ ). In the same hypothesis, 
we find a median value of agesF/V = 0.1 for the FULL 
sample and a similarly low value for all the sub-samples. 
However, the large uncertainties on agesF and r reflect 
in a wide range of allow ed agesF It ratios for all SFHs, as 
shown in Tab. [7]andj8 Erb et al. (2006) used the expo- 
nentially decreasing SFR model to find the best fit SED 
of individual star forming galaxies at z ~ 2. As they 
found that the r parameter space was basically allowed 
in its entirety, they chose the star formation history that 
allowed them to recover SFR values in agreement with 
Ha SFR estimations. We are now motivated to estimate 
the SFR averaged over the last 100 Myr (Equation|8| and 
compare it with SFR corr (UV) and SFR. We will present 
it and comment about the (dis) agreement in the next 
section. 

5.1.3. Star formation rate estimations 

The instantaneous SFR is very sensitive to values of 
both agesF and r. We therefore expect high uncertain- 
ties in its estimation as reported in Tab. [5] - [HJ In the 
majority of the cases we report its value for complete- 
ness, even if the large error bars do not make it a use- 
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ful quantity to compare with the literature. Assuming 
a CSF, we expect the instantaneous SFR, the averaged 
SFR over the liletime ol the galaxy, and the SFR corr to 
agree. However, if a galaxy is less than 100 Myr old, 
<SFR> 10 o would be lower than the SFR averaged over 
the lifetime of the galaxy. The latter is expected to be 
smaller than the instantaneous SFR by a factor roughly 
corresponding to 100 Myr/agesF- This behavior can be 
observed for instance in the FULLsample for each of the 
three SFHs (Tab. [5]), where the instantaneous SFR is 
higher than <SFR> 10 o within the 68% confidence level 
range. On the other hand, SFR corr (UV) = 15 [6-30] M 
yr -1 is roughly consistent with <SFR>i 00 = 4 [2-9] M© 

yr _1 _ 

A similar behavior is observed in the sub-samples. For 
the ones with younger SEDs (IRAC bright, red-LAEs, 
UV-bright), <SFR> 100 and SFR corT .(UV) imply lower 
values of star formation rate than the instantaneous SFR. 
It is interesting to note that for red-LAEs, the inferred 
values of instantaneous SFRs are as high as thousands 
M yr _1 . For the sub-samples with older SEDs (IRAC- 
faint, UV-faint, high-EW, where the agesF can reach sev- 
eral hundreds of Myr) , the three SFR indicators are con- 
sistent 

In j j5.4| we will compare the SFR estimations of the 
FULLsample with that inferred by X-ray emission. 

5.2. LAE colors 

The red (B — R> 0.5) sub-sample contains only 15% 
of our sample LAEs, indicating that the majority of our 
z ~ 2.1 LAEs have blue rest-frame colors. The red-LAE 
sub-sample has the highest IR flux of any sub-sample as 
we can see in Fig. [4] A more global description of the 
UV-through-IR SED is offered by the R- [3.6] (rest-frame 
NUV- J) color which is sensitive to the UV slope and the 
Balmer/4000A break, so it traces a combination of dust 
reddening and age. The 3.6 /xm magnitude is a tracer of 
the stellar mass. We plot the R- [3.6] color versus the 
3.6 fj,m band magnitude to study the color— stellar mass 
trend of the z ~ 2.1 LAEs and their sub-samples in Fig. 
[9^,. It appears to be a linear relation between the mass 
and this color, with the IRAC-faint, UV-faint, and high- 
EW also being the bluest and least massive sub-samples. 
The red-LAEs appear at the upper left corner of the fig- 
ure being a special sub-sample of extremely massive and 
dusty galaxies. They can be even dustier than spectro- 
scopically confirmed z ~ 2.1 BX galaxies although their 
stellar masses can be similar. The BX are brighter in R 
due to photometric and spectroscopic selection, so they 
lie outside the linear mass-color relation. L08 found a 
continuum of properties in this color-magnitude plot be- 
tween z ~ 3 LAEs and LBGs, observing that z ~ 3.1 
LAEs were at the faint, blue end. We can say that 
z ~ 2.1 LAEs have a wide range of masses and there 
is a continuum of properties from the redder and more 
massive to the bluer and less massive sub-samples. Also, 
the datum corresponding to L08und sub-sample is lo- 
cated outside the linear mass-color relation, showing the 
less dusty and/or younger nature of the typical z ~ 3.1 
LAEs than those at z ~ 2.1. 

In Fig. [9J) we show the R- [3.6] color versus the R 
magnitude, which is roughly proportional to the SFR 
derived from the rest-frame UV continuum assuming no 



dust extinction. The intrinsic SFRs calculated at z ~ 2.1 
are much higher than at z ~ 3.1 but z ~ 2.1 LAEs 
are dustier. In this figure the dustiest red-LAEs are as 
faint as the IRAC-bright and UV-bright sub-samples in 
the rest-frame UV continuum. Therefore, we see a less 
clear linear mass-SFR relation. The BX galaxies are the 
brightest objects as, by definition, they are selected to 
have R < 25.5. From this figure we can also see that 
the R band sensitive SFR is comparable in the L08det, 
L08und, and our z ~ 2.1 LAE sub-samples, but the UV- 
bright, IRAC-bright, and red-LAEs are characterized by 
higher SFR values. 

5.3. Dust properties in LAEs 

In this sub-section, we investigate the dependence 
of the dust on photometric properties, such as l^Lyu, 
EW(Lya). 

We find that wherever the dust amount is signif- 
icant, the Lya EW tends to be low. In fact the 
median(EW(Lya)) = 30 A for the typical IRAC-bright 
(E(B-V) = [0.3-0.5]), while it is 70 A for the typical IRAC- 
faint ((E(B-V)=[0.0-0.3])) galaxy. Galaxies brighter in 
the continuum, with M* > 10 9 M and characterized by 
a significant amount of dust, show low equivalent width 
Lya lines. This result seems to be in disagreement with 
the hypothesis of enhancement of EW(L yo) due to the 
presence of a clumpy in terstellar medium ( Neufeld||l991 



Finkelstein et al . (2 008 ) . In this scenario, dust grains en 



closed inside clumps of HI would not be able to absorb 
Lya photons, as the HI density gradient would scatter 
them away from clouds, favoring their escape from the 
galaxy. 

We relate the dust amount to the rest-frame UV slope 
and we also compare these quantities with the equiva- 
lent width and Lya luminosity. The slope of the rest- 
frame UV spectrum is proportional to the reddening of 
the stellar continuum measured by the E(B-V) param- 
eteip*] we approximate the flux density as a power law 
(t\ tx \^ uv ) and estimate its index, Puv- We fit the 
rest-frame UV spectrum from B to I to obtain f}jjy 
and its uncertainty. We use Calzetti's assumption that 
A(1600 A) = 4.39 x E„(B-V), E(B-V) = (0.44±0.03) 
x Eg(B-V), and the |Meurer et al.||1999| relation A(1600 
A) = 4.43 + 1.99 x p vv . Therefore, high-E(B-V) ob- 
jects are characterized by high, even positive, Pjjv val- 
ues. E(B-V)=0 corresponds to f\ cx A -2 2 , nearly flat in 

In Fig. [TO] we show the fiuv index as a function of 
Llj, q and EW res t-/rame- From the left panel we can see 
that for higher values of Lya luminosity {Lhya > 10 42 
erg sec -1 ) j3uv tends to -1.9 with a scatter of 0.2. This 
means that high- luminosity LAEs at z ~ 2.1 present a 
flat rest-frame UV spectrum. In the low-luminosity bin 
$xjy=-\.h and the scatter is much larger (p-p — 1.1). This 
is due to the low L L ya objects that can be characterized 
by Puv up to 3. Tracing a horizontal line at fiuv = - 
1.9 and a vertical line at log(Lij, Q ) = 42.35, 62% of the 
less luminous and 45% of the more luminous objects have 
Puv > -1-9. The right panel shows that the highest-/3[/v 

15 Note that our quantity E(B-V) is referred to as E S (B-V) by 
Calzetti et al. (2000), while the color excess derived from the 
nebular gas emission lines is defined as E 9 (B-V). 
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sources are red-LAEs, with EW res t_/ rame < 50 A and 
steeper spectra, mean(/3^y) = 0.17±1.15. Also, the UV- 
bright LAEs present EW res t_/ rame < 100 A and j3uv = 
— 1.2±0.95, while the objects brightest in Lya luminosity 
are characterized by a flatter spectrum (/3jjv = — 1-9 ± 
0.3) with small scatter in /3jjv and EW„ S [_/ rame up to 
200 A. In the same figure we also show the data from 
Nilsson et al. (2009) at z ~ 2.3. In their sample, LAEs 
with the highest EW values are also characterized by 
flat rest-frame UV spectra and high Lya luminosities. 
Steeper spectra are seen for z ~ 2.3 EW reS i_/ rame < 40 
A objects, wh ich have lower Lya luminosity. 

In Fig. [TT] we plot the R magnitude versus the stel- 
lar continuum dust reddening, E(B-V), estimated from 
the ratio between the observed SFR(UV) and SFR(Lya) 
from GulO for the individual LAEs of our z ~ 2.1 sample. 
The observed SFR(UV) and SFR(Lya) are related to the 
intrinsic SFR by assuming the dust law from Calzetti 
et al. (2000) for stellar continuum and nebular emis- 
sion. We assume that the stellar continuum reddening is 
proportional to the nebular reddening, (E(B-V) = c x 
E g (B-V)), with the c value depending on the dust ge- 
ometry of the ISM i.e., dust enhancement surrounding 
star forming regions. The rough linearity in the plot oc- 
curs because SFR(UV) is obtained from the observed 
R band flux density. We first assume Calzetti's con- 
stant c = 0.44, obtained from the study of local star- 
burst galaxies. The median(E(B-V)) is 0.06 for the full 
sample. Some recent results suggest that dust properties 
at higher redshifts are differ ent and a c ~ 1 is required 
to explain the observations (Erb et al. 2006a). There- 
fore, we also plot R versus reddening assuming c = 1.0. 
In this case the median(E(B-V)) is 0.3 for the full sam- 
ple. Comparing with the E(B-V) inferred by our SED fit 
(E(B-V) FUL Lsa m pio = 0.22[0.00-0.31], E(B-V) UV -w = 
0.04[0.00-0.24], E(B-V) UV -bright = 0.32[0.09-0.38]), pre- 
sented in the plot as black triangles, we investigate which 
is the proper factor of proportionality between E(B-V) 
and Eg(B-V) in the framework of the Calzetti law. Un- 
der the assumption of simple radiative transfer of Lya 
photons, Fig. 11 is compatible with any factor of pro- 
portionality 0.44 < c < 1.0 at z ~ 2.1, but c = 1.0 is 
disfavored at 68% confidence level. Uncertainties in Lya 
radiative transfer make this comparison even less strin- 
gent, but a future independent determination of c would 
allow this comparison to be used to constrain radiative 
transfer effects. 

5.4. X-ray stacking constraints on SFR and obscured 

AGN 

In the central part of our field, the CDFS, we use the 
4 Msec Chandra data to estimate the average stacked X- 
ray flux among 38 FULLsample objects. We obtain an 
upper limit in the Chandra soft-band flux and a > 2er de- 
tection in the hard band. The rest-frame hard X-ray flux 
is sensitive to the presence of high-mass X-ray binaries 
(HMXBs) in the galaxies. In this case In this case the 
involved time scale is 10-100 Myr, in which O and B stars 
compose the binary systems. For longer periods of time, 
low-mass X-ray binaries (LMXBs) and thermal emission 
can also contribute to the total hard X-ray b udget, even if 
at z ~ 2 the number of LMXBs is negligible ( Ghosh et al 



X-ray data as anoth er way to infer the SFR (see also Kur 



2001 1. The HMXBs trace star formation, so we use these 



czynskiet al. 20101. The 4 Msec data are so deep that 
only very- low luminosity or heavily-obscured AGN would 
escape individual detection. As the observed-frame soft 
band could be obscured, we concentrate on the observed- 
frame hard ba nd to infer the SFR . 

We use the IPersic et al.M2004l indicator to turn the 
hard X-ray luminosity into SFRp ers ; c . They assume that 
/ = (HMXB luminosity) /(total hard X-ray Luminosity) 
equal to 1 for high redshift systems implies that the hard 
X-ray l uminosity is entirely from the HMXBs. We also 
use the |Lehmer et al.|[20l0] equation that takes into ac- 
count the contribution of LMXBs together with that of 
HMXBs to estimate the SFR, h X -ra y = L(LMXB) + 
L(HMXB) = a M* + /3 SFR). This way they fit the 
data with less scatter than Persic across a wider range 
of SFRs and X-ray luminosities. When neglecting the 
factor involving the low-mass X-ray binary luminosity, 
the Lehmer and Persic equations are consistent. Our 
LAEs are characterized by low values of M* , making the 
"a M*" term negligible with respect to the 11 /3 SFR" 
term. In Tab. [9] we present stacked observed-frame soft 
and hard X-ray fluxes, and the inferred SFRs. The er- 
rors are estimated by measuring the background in ~ 
30 apertures around the central region, with the same 
size as for the source. Then, we calculate the mean and 
standard deviations of those background measurements. 
Obtaining an upper limit means that the stacked signal 
is not detected at the 2a level. The observed-frame soft 
band X-ray flux is < 2.2E-18 erg sec -1 cm~ 2 , while the 
observed-frame hard band X-ray flux is (4.6 ± 1.9)E-17 
erg sec" 1 cm" 2 . The inferred SFR Perslc is 1180 ± 480 
M yr" 1 , and SFR Lehmer is 730 ± 300 M yr" 1 , by as- 
suming a = and (3 = (1.62 ± 0.22) E+39 erg sec -1 
(Mo/yr)" 1 . 

Our SED fitting analysis shows that the typical LAE 
has an instantaneous SFR of 35 [0. 003-170] M yr" 1 and 
the averaged <SFR> 100 is 4 [2-30] M yr _1 . Over the 
last 100 Myr and for a CSF assumption, SFR is directly 
comparable to the rest-frame UV flux density. Tak- 
ing into account the dust amount through the Calzetti 
et al. (2000) law, we derive SFR corr (UV). Based on 
the E(B-V) value measured in the SED fitting using a 
CSF, SFR corr (UV) turns out to be 15 [6-30] M yr" 1 
and <SFR> 100 is 4[2-9] M yr" 1 . Hence, for the 
FULLsample, SFRs of ~10 M yr -1 seem to be de- 
rived by <SFR> 10 o, SFR corr (UV), and instantaneous 
SFR estimations, independently of the SFH, but the 
SFR /lard x-ray is significantly higher. 

There are a few plausible explanations for this 
discrepancy. We could be seeing a diffuse thermal 
plasma in X-rays. Persic fc Rephaeli] (2002) discuss the 
diffuse component of X-ray luminosity in the context 
of X-ray SFR calibration, specifically. The discrepancy 
could also be due to a contamination by a small 
number of low- luminosity or heavily obscured AGNs 
in our sample. They would dominate the averaged 
X-ray stack, but the median approach we are using 
in the SED analysis is not sensitive to them. We 
look at the hardness ratio (HR=(countsxray-/tard- 

COUntS X ray-soft) / {cOUntS X ray-hard+COUIltSXray-soft)) 

and at the ratio between X-ray and R band fluxes to try 
to identify a discriminator between AGNs and pure star- 
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forming galaxy imprints in our X-ray stack. Looking at 
Fig. 5 of Treister et al. (2009) and Fig. 8 of Wang at al. 
(2004), HR > 0.43 and \og{¥ xray-hard/? r) = -0.8±0.4 
values are consistent with having a mixed population of 
obscured AGNs and pure star-forming galaxies in our 
LAE sample. 

6. CONCLUSIONS 

In this work we calculated the observed spectral en- 
ergy distribution of the z ~ 2.1 LAE sample and its 
sub-samples (Tab. |2j), obtained based upon photometric 
properties by dividing the sample roughly in halves. We 
mainly concentrated on the IRAC-bright and faint sub- 
samples, because the separation at observed 3.6 /^m is 
expected to be a separation based on the stellar mass; 
on the red-LAEs as a special category with observed 
SED different from the bulk of z ~ 2.1 LAEs; the UV- 
bright and UV-faint sub-samples, because the separat- 
ing R magnitude represents the rest-frame UV contin- 
uum of our galaxies. The complete z ~ 2.1 LAE sam- 
ple, defined in GulO and here in §2, is used to investi- 
gate the EW correlation with SED parameters. To al- 
low the properties of our LAE sample to represent the 
general properties of LAEs at this redshift in the Uni- 
verse, we take sample variance into account in our error 
estimation by using bootstrapping technique (§2). We 
used the CB10 stellar population code to generate con- 
tinuum SED models, probing a wide range of star forma- 
tion histories parametrized by the e-folding time, r. We 
assumed the Salpeter IMF to be able to compare with 
results in the literature. It was also justified by the fact 
that the stacked flux density in the NB5015 produced a 
non-detection (<f3l). In fact at z ~ 2.1 that narrow-band 
flux would reveaTthe Hell emission, imprint of top heavy 
IMF and/or very-low metallicity stellar population. The 
EW(HeII A1640), implied by NB5015 non detection, is 
consistent with zero, as it may be expected for a normal 
IMF galaxy. We looked for the best fit SED (Fig. |J Fig. 
|4) and the 68% confidence level of fitting parameters of 
the typical LAEs of our sample and sub-samples. In Fig. 
[3 [5 [6j and [7] we presented joint contours of different 
combinations of parameters and star formation histories. 
Equally good fits of our data can be obtained assuming- 
exponential increasing, decreasing or constant star for- 
mation rates. The stellar mass, M* , and the reddening, 
E(B-V), are well constrained, while agesF has large un- 
certainties and r is not constrained at all. 

Heterogeneity can be revealed in a sample of galaxies 
either by fitting ind ividual galaxies, if ther e is enough 



signal-to-noise (e.g., Finkelstein et al. 2008), or by our 



approach of stacking sub-samples selected according to 
their observed photometric properties. Our results are 
in agreement with Nilsson et al. (2009, 2011). The GulO 
clustering result showed that z ~ 2.1 LAEs evolve into lo- 
cal Universe galaxies with median luminosity equal to L* 
and that z ~ 2.1 LAEs therefore represent the building 
blocks of z ~ Milky Way-type galaxies. From our me- 
dian statistics analysis, our results show that these LAEs 
appear characterized by median M* < 10 9 M© and a ~ 
40 Mq yr _1 median instantaneous SFR. This moderate- 
to-high SFR value might be caused by rapid accretion 
of mass, which, together with merging phenomena, must 
play a major role in causing LAEs to evolve into z ~ 
L* galaxies. The observed starburst is presumably not 



the only phase of active star formation happening in- 
side these galaxies between formation and the present 
day; multiple episodes of SF are likely to contribute to 
the growth of stellar mass. We could also be observing 
Lya emission as the effect of the SF in particular loca- 
tions within a galaxy where the Lya photons are free to 
propagate, even if other galaxy regions are temporarily 
inactive or dustier. 

The main results obtained from SED fitting and fur- 
ther discussions are the following: 

• z ~ 2.1 LAEs appear to be low mass galaxies with 
moderate amounts of dust. The instantaneous SFR val- 
ues are estimated of the order of tens of M Q yr _1 , with 
an extreme value of 170 M© yr _1 within the 68% confi- 
dence level. However, the star formation rate averaged 
over the last 100 Myr, <SFR>ioo, is consistent with ~ 
4 M Q yr -1 . The mean stellar age, t*, is ~ 10 Myr, so 
that the LAE SED is consistent with being dominated by 
young stars i.e., LAEs are undergoing a period of intense 
star formation. 

• z ~ 2.1 LAEs appear in median dustier (E(B-V) = 
0.22[0.09-0.28]) than z ~ 3.1 LAEs assuming constant 
star formation rate model, but present consistent stellar 
masses. 

• The typical z ~ 2.1 LAEs appear to be less massive 
(M* < 10 8 M ) than z ~ 2.1 UV-continuum selected 
star forming galaxies (Fig. [8]), implying that LAEs tend 
to occupy the low-mass ena of the distribution of star 
forming galaxies at z ~ 2. 

• The IRAC-bright LAEs are in median significantly 
more massive and dustier than those that are faint in 
IRAC. The R < 25.5 LAEs are also more massive than 
the fainter ones. Our typical z ~ 2.1 LAE is blue with 
median(_B — R) = 0.13, but the 15% of the objects clas- 
sified as red-LAEs are significantly more massive and 
dustier, with E(B-V) > 0.3. 

• The sub-sample of LAEs characterized by high- 
EW(Lya) shows the lowest stellar mass and a moderate 
dust amount (E(B-V) = [0.0-0.3]). 

• z ~ 2.1 LAEs have a wide range of stellar masses 
(Fig. [9]). There is a continuum of properties from the 
redder and more massive to the bluer and less massive 
typical LAEs, which are the majority. This is another 
indication that LAEs are located at the lower mass end 
of the 2^2 star forming galaxy mass function, where 
the bulk of galaxies lie. 

• The red-LAEs are characterized by EW res t-/rame < 
50 A (FigfTol). The Lya luminous LAEs are consistent 
with flat rest-frame UV spectra and can present equiva- 
lent widths up to 200 A. 

• The best-fit instantaneous SFR and <SFR>ioo differ 
significantly in several sub-samples, making it important 
to consider the timescales of the various SFR estimators 
when studying starburst galaxies, such as LAEs. 

Our results point out differences in the dust amount 
and SFR of LAEs between redshifts z ~ 3.1 and z ~ 2.1. 
A possible explanation is that the redder and dustier 
LAEs, observed at the later time in the evolution of the 
Universe, are experiencing evolved phases of star forma- 
tion, unhke_WTfih_dus^free z ~ 3.1 counterparts (Ga07, 
L08, Acquaviva et al.|2011 ). Future work will investigate 
whether these results are robust by using uniform SED 
fit methodology and stellar population synthesis model- 
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ing for the two epochs. If the indications coming from 
this work and N10 are confirmed, this will be a key step 
in our understanding of how the building blocks of L* 
galaxies evolve with cosmic time. 
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TABLE 1 

Limiting magnitude in the 13 MUSYC broad bands in ECDF-S 



u 


B 


V 


R I z 


J 


H 


K 


3.6/im 


4.5/im 


5.8/im 


8.0^m 


26.1 


26.9 


26.5 


26.5 24.6 23.6 


22.8 


22.3 


21.5 


23.9 


23.7 


21.9 


21.8 



TABLE 2 
Sub-sample definitions 



I Sub-sample name Selection Nlae 



FULLsample 


entire z ~ 2.1 LAE sample 


216 


IRAC-bright 
IRAC-faint 


hstim > 0.57/^J 
hsfj-m < 0.57/^J 


47 
67 


red-LAE 
blue-LAE 


B - R > 0.5 
B - R < 0.5 


34 
182 


UV-bright 
UV-faint 


R < 25.5 
R > 25.5 


118 

98 


high-EW 
low-EW 


log(L iyce ) > 42.1 
EW reat _/ rame >66 A 

log(L £va ) > 42.1 _ 
EW rest _/ rame <66 A 


60 
59 



Note. — The number of LAEs in each sub-sample is based on 
having coverage in the optical bands. The only exception is for 
the sub-samples defined based on IRAC 3.6 /J.m, which arc selected 
from the 114 LAEs in the IRAC "clean" regions. 



TABLE 3 

Overlap of objects among the sub-samples 



sub-sample 


IRAC-bright 


IRAC-faint 


red-LAEs 


blue-LAEs 


UV-bright 


UV-faint 


EW-high 


EW-low 


IRAC-bright 


47 





16 


31 


42 


5 


4 


14 


IRAC-faint 





67 


5 


62 


16 


51 


30 


14 


red-LAEs 


16 


5 


34 





25 


9 


8 


6 


blue-LAEs 


31 


62 





182 


93 


89 


52 


53 


UV-bright 


42 


16 


25 


93 


118 





5 


49 


UV-faint 


5 


51 


9 


89 





98 


55 


10 


EW-high 


4 


30 


8 


52 


5 


55 


60 





EW-low 


14 


14 


6 


53 


49 


10 





59 



NOTE. — For each sub-sample wc report the number of objects belonging to other sub-samples. Boldface numbers indicate the actual number 
of objects in each sub-sample. 



TABLE 4 

Free parameters and value ranges used in the SED 
fitting process 



parameter 


Min. Value 


Max. Value 


step 


Parameter 


Min. Value 


Max. Value 


step 


E(B-V) (mags) 


0.00 


0.60 


0.01 


log(M,/M ) 


6.00 


12.00 


0.05 


age SF (Gyr) 


0.005 


3.00 


~0.1 dcx 


r(Cyr) 


-4 


18 


~0.15 dex 



TABLE 5 
Full sample SED parameters 



SFH 


log(M»/M ) 


E(B-V) 


agesF(Gyr) 


r(Gyr) 


SFR (M /yr) 


<SFR> 100 (M /yr) 


Xbest / lt / d -°' f 


all 


8.6[8.4-9.1] 


0.22[0. 00-0.31] 


0.018[0.009-3] 


-0.02 [all] 


35[0.003-170] 


4 [2-30] 


12.3/9 


r > 


8.5[8.4-9.2] 


0.20[0. 00-0.29] 


0.010[0.009-0.227] 


0.007[all] 


14 [0.003-60] 


3[2-15] 


12.4/9 


r < 


8.6[8.5-9.0] 


0.22[0.07-0.31] 


0.018[0.009-3] 


-0.020[all] 


35 [5- 170] 


4 [2-30] 


12.3/9 


CSF 


8.6[8.4-9.0] 


0.22[0.09-0.28] 


0.012[0.009-0.161] 


oo 


30 [6-60] 


4 [2-9] 


12.6/10 



NOTE. — M* is the observed stellar mass, taking into account the gas reprocessing inside a galaxy. agesF is the age since star formation began. 
SFR is derived from the star formation history as y^j-exp (- ag l SF ) (Equation 13. For the constant star formation model (CSF) wc fix r to be an 
arbitrarily high value of 18 Gyr in our code. 
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TABLE 6 

2 ~ 2 AND Z ~ 3 SED FITTING 



reference 


z 




E(B-V) 




a g e SF(Gyr) 


OPD (~\A^ /\rr\ 


N07 


3.15 


oo 


0.07[0.02-0.09 


8.67[8. 18-8.95] 


0.85[0.43-0.98] 


0.66[0.35-1.17] 


Ga07 


3.1 


0.75±° f 5 


0.00[0.00-0.08 


9.00[8. 80-9.20] 


0.02[0.01-0.05] 


2.00[1.00-3.00] 


OlOund 


3.1 


oo 


0.03[0.00-0.08 


8.12[7.75-8.38] 


0.07[0.01-0.18] 


0.36[0.19-0.78] 


OlOdet 


3.1 


oo 


[0.00-0.70] 


[8.97-10.43] 


[0.005-0.407] 


[11-5600] 


N10 


2.3 


2 SSP 


0.07 


0.03-0.20 


9.83[9. 66-9.98] 


0.07[0.02-0.09] 




FULLsample (CSF) 


2.1 


oo 


0.22 


0.09-0.28 


8.6[8.4-9.0] 


0.012[0.009-0.161] 


30 [6-60] 


FULLsample (r > 0) 


2.1 


0.007 [all] 


0.20 


0.00-0.29 


8.5[8.4-9.2] 


0.010[0.009-0.227] 


14[0. 003-60] 


FULLsample (t < 0) 


2.1 


-0.020 [all] 


0.22 


0.07-0.31 


8.6[8.5-9.0] 


0.018[0. 009-3] 


35[5-170] 


IRAC-bright 


2.1 


oo 


0.40 


0.32-0.46 


9.4[9.2-9.5J 


0.009[0.006-0.015] 


260[110-530] 


L08dct 


3.1 


oo 


0.00 


0.00-0.10 


9. 95[9. 78-10.08] 


1.60[1.20-2.00] 


6 [5-7] 


IRAC-faint 


2.1 


oo 


0.11 


0.00-0.25 


8.3[8.0-8.9J 


0.045[0.009-0.716] 


5[l-25] 


LOSund 


3.1 


oo 


0.00 


0.00-0.10 


8.48[8. 00-8.85] 


0.16[0.05-0.30] 


2 [1-3] 


BX 


2.0 %spec 2.2 


oo 


0.35 


0.31-0.37 


9.6[9.5-9.7J 


0.009[0.009-0.01] 


470[315-526] 



Note. — The references indicate Nilsson ct al. (2007)=N07, Nilsson ct al. (2011)=N10, |Ono et al.| (2010)=O10 with K(3<r)=24.1, 
Gawiscr et al. (2007) — Ga07 with mag^e [3-6] > 25.2 and a 2-populat ion model. CSF indicates constant star formation. L08det, L08und 
arc IRAC-detected and IRAC-undetected subsamples from [Lai et al.] (2008). FULLsample, IRAC-bright, IRAC-faint, and BX arc from 
this work. 
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TABLE 9 

X-RAY STACKING RESULTS FOR THE FULLSAMPLE 



obs soft 
fluxx _ ray 


obs hard 

fluiix-ray 


rest hard 


rest hard 

— ray 


SFR, Persic 


SFR,£ e /,, meT . 


SFR corr (UV) 


<SFR>i 00 


SFR 


<2.2E-18 


(4.6±1.9)E-17 


(3.8±1.5)E-17 


(12±5)E+41 


1180±480 


730±300 


15[6-30J 


4[2-30J 


35[0.003-170] 



Note. — The units are erg s 1 cm 2 for the observed-frame soft-band (obs soft fluxx — ray), observed-frame hard-band (obs hard fluxx-ray). and 
rest-frame hard-band (obs hard fluxx — ray) X-ray flux, erg s~ 1 for the rest-frame hard-band X-ray luminosity (rest hard Lx_ ray ), and Mq yr _1 for 
the star formation rate estimates. The dust-corrected SFR(UV) corresponds to the values calculated assuming a CSF model, while the instantaneous 
and averaged SFR considering the entire parameter space. 




log(« (A) 



Fig. 1.— Filter transmission curves for MUSYC ECDF-S broadband photometry. From left to right: U, B,V, R, I , z, J, H, K , IRAC 3.6, 
4.5, 5.8, 8.0Atm. 




Fig. 2. — Stacked SED of the FULLsample of z ~ 2.1 LAEs (black triangles with error bars). In solid lines we show the spectral energy 
distribution corresponding to the exponentially decreasing (magenta), exponentially increasing (red), and constant (orange) star formation 
histories, generated using the best fit parameters reported inside the plot. Dashed lines illustrate four models at the edge of the 68% 
confidence level. The fit parameters for these models are M„=8.91 X 10 8 M Q , E(B-V)=0.07, age SF =0.2273 Gyr, t=-1.39 Gyr (cyan); 
M*=10 9 M , E(B-V)=0.1, age SP =0.114 Gyr, t=0.118 Gyr (brown); M,=2.82 X 10 s M Q , E(B-V)=0.16, age SF =0.018 Gyr, r=-4.0 Gyr 
(blue); M„=2.82 X 10 s Mq, E(B-V)=0.17, age SF =0.018 Gyr, r=2.812 Gyr (green). 
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Fig. 3. — 68% confidence levels for the model parameters in the case of the exponentially decreasing (blue), exponentially increasing (red), 
and constant (green) SFHs. Different panels show different combinations of parameters. The contours show that while M* and E(B-V) 
are well constrained, a wide range of ages and even wider range of t are allowed at 68% confidence. Fits of comparable quality can be 
obtained for all three SFHs. The best fit models are indicated as colored dots. The dot with the black circle represents the best fit in the 
entire parameter space. 
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Fig. 4. — Observed median SEDs (color symbols) and best fit model SEDs (color curves) in a log-log plot. The FULLsample best fit SED 
is plotted (black dashed curve) for comparison. Each panel shows pairs of sub-samples as labelled. 
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Fig. 5. — 68% confidence level in Log(M*/M0) versus E(B-V) plane for the exponentially decreasing (blue), exponentially increasing 
(red), and constant (green) SFHs for different LAE sub-samples. The dots represent the best fits for each of the three SFHs. The dot 
with the black circle represents the best fit in the entire parameter space. These parameters are well constrained for every combination of 
sub-samples and SFHs. More massive sub-samples are characterized by higher dust reddening. 
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Fig. 6. — 68% confidence level in E(B-V) versus log(agegp) plane for the three SFHs, with the symbols and colors as in Fig. \5\ We sec 
minimal degeneracy between these two parameters. agegF is not well constrained except that IRAC-bright and UV-bright suD-samples 
have ages less than 20 Myr in the exponentially declining and constant SFHs. 
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Fig. 7. — 68% confidence level in Log(|r|) versus log(agegp) plane for the three SFHs, with symbols and colors as in Fig. \5\ The 
contours show that even if agesF and r are not well constrained individually, some regions of the parameter space are disfavorea at the 
68% confidence level. 
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Fig. 8. — Observed SED (brown circles with error bars) for the sample of spectroscopically confirmed BX star forming galaxies at 
2.0 < z spec < 2.2. The solid brown curve shows our best fit SED model, obtained with an exponentially decreasing SFH, while the solid 
blue and dotted cyan curves show the best fit obtained assuming constant and exponen tially decre asing SFR. The dashed magenta line 
shows the SED which corresponds to the median best fit parameters of the sample from Erb et al. 2006b As a comparison we show the 
best fit model SED of our UV-bright sub-sample in olive green. 
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Fig. 9. — Color-magnitude diagram for color-mass (left) and color-SFR (right) relations. The symbols represent the z ~ 2.1 LAE sub- 
samples with the color coding as in Fig. W] (see legend). The black stars represent the results of Lai et al. (2008) for the z ~ 3.1 IRAC 
L08und (on the extreme right) and L08det LAE sub-samples. The larger horizontal error bars in the left panel are caused by larger flux 
uncertainties in the observed IR in comparison to the optical bands. 
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Fig. 10. — Rest-frame UV slope (Puv) versus hhya (left) and rest-frame equivalent width (right) for z ~ 2.1 LAE samples. The y-axes 
also show the stellar continuum E(B-V) parameter from the |Meurer e t al. 1999 formalism. Black dots with gray error bars correspond to 
the z ~ 2.1 LAEs. In the left panel we show the z ~ 2.3 data separating them based upon equivalent width: EW re3t _ f rame < 40 A in 
red, 40 A < EW rest _ frame < 80 A in orange, and EW rest _/ rame > 80 A in blue crosses. In the right panel the orange crosses with error 
bars correspond to the z ~ 2.3 LAEs from Nilsson et al. (2009). Three z ~ 2.1 sub-samples are also presented here: red LAEs as red dots, 
LAEs with l^Lya > 10 42 ' 3 erg sec -1 as magenta circles, and the UV-bright LAEs as small green dots. 
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Fig. 11. — Observed R band magnitude versus stellar continuum reddening, E(B-V), inferred from the Lya/UV luminosity ratio and the 
assumption of simple Lycr radiative transfer. We assume the dust law from Calzetti et al. (2000) and E(B-V) = c X E g (B-V). Red dots 
correspond to c = 0.44 and blue dots to c = 1.0. Black triangles represent the median observed R magnitude and the best fit E(B-V) values 
from the SED analysis for the UV-bright and UV-faint sub-samples. 



